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Abstract 



o 
o 

(N 

The spectra of relic gravitational waves produced as a result of cosmological expan- 
sion of the inflationary models are derived in Brans-Dicke (BD)theory of gravity.The 
time dependence of the very early Hubble parameter and matter energy density are 
derived from frequency dependent spectrum of relic gravitational waves. Also it is found 
that Brans-Dicke scalar field contributes to the energy density of relic gravitons. 
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^O 1 INTRODUCTION 

O 

Relic gravitational waves are, quite possibly, the only source of information on the physical 
conditions in the very early universe. The energy density and spectrum of relic gravitational 
waves depend on a specific rate of expansion of the universe in its very distant past. The ac- 
Vh ! tual detection of relic gravitons, or even the experimental restriction of their possible spect 

ral energy density, are capable of producing fairly meaningfull conclusions about the param- 
eters of the very early universe [1]. In general, the information stored in the relic gravitons 
allows one to judge the overall evolution of the cosmological scale factor. However, in more 
specific models of the universe, this information gives us the direct data on the variability 
of the Hubble parameter. From the theoretical point of view, the most natural alternatives 
to General relativity(GR) are the scalar-tensor theories, which contain an additional scalar 
field, 0, whose relative importance is determined by the arbitrary coupling function uj{4>) 
[2].GR corresponds to the case where the field 0=const.=l implying an infinite uj((f)).The 
Brans-Dicke theory is defined by a constant oj ((f)) [3] .This is a successful theory because 
it explains almost all the important features of the evolution of the universe. Some of the 
problems like inflation, early and late time behaviour of universe [4], cosmic acceleration and 
structure formation [5], cosmic acceleration, and coincidence problem[6] can be explained 
in the B-D formalism. For large u;,B-D theory gives the correct amount of inflation and 
early and late time behaviour, and for small negative u it correctly explains cosmic accel- 
eration, structure formation and coincidence problem. Also, this theory appears naturally in 
supergravity theory, Kaluza-Klein theories and in all the known effective string actions and 
thus, in itself, justifies further consideration. 

The fundamental feature of BD theory is that there is a scalar field (j) coupled to the 
curvature. Its essential role in homogeneous cosmologies is to induce a time variation of the 
gravitational constant. Also this scalar field <p modifies the expansion rate in the early uni- 
verse. Cosequently, flatness problem, horizon problem, the primodial nucleosynthesis prob- 
lem, and the energy density and spectrum of relic gravitational waves are all modified. Here 
we aim to assess the modifications induced by the coupling of the scalar field to the en- 
ergy density and spectrum of the graviton. Previous work on the subject of gravitational 



radiation within Brans-Dicke theory, carried out by Wagoner [7] was concerned with gravi- 
tational waves generated in the weak field limit. Later, Barrow[8] derived the gravitational 
wave equation for the general class of scalar tensor theories and solved for vacuum and 
radiation era for different cosmological models. Also Sahoo and Singh[9] solved the same 
wave equation for all epochs of the universe and found that for negative present values of 
the parameter to gravitational waves grow in all epochs of the universe. 

An outline of the paper is as follows. In the next section we write the solutions of the 
gravitational wave equation in BD theory .The Bogoliubov transformation coefficients and 
energy density of gravitons are derived in the section three. In the section four frquency 
dependence of early Hubble parameter and matter energy density are derived. We conclude 
with section five. 



2 BASIC EQUATIONS 

The perturbations hij to the metric that represent weak gravitational waves can be expressed 
as 

h(t,x) = d nh n (t,x) 



h n (t,x) 



a 2 (t) 



Y n (t)U n (x) 
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The wave number ,n,n 2 = (rii) 2 + (ri2) 2 + (n^) 2 ,is associated with the frequency ^(measured 
in Hz) according to the relation v = -£— in units c = G = % = l.The functions U n and Y n 
satisfy 
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In the above equations a(t) is the scale factor, V 2 is the spatial Laplacian operator, an 
overdot indictes derivative with respect to the synchronous cosmic time t.Equation(4) is 
gravitational wave equation for zero curvature FRW Universe [8]. The Brans-Dicke scalar- 
tensor theory of gravity is characterised by the restriction that uj(4>) = lu = constant and 
U = dU /d<p = O.We introduce the conformal time r] defind by 

a(r])dr] = dt. (5) 

Now each component h n of the gravitational wave perturbations can be written as 



h n (r],x) 



1 



R(v) 



^n(v)Un(x), 



(6) 



where, 

R( V ) = a( V )4> 1/2 ( V )- (7) 

So, for the Brans-Dicke scalar-tensor theory of gravity and for Eq(6), the gravitational wave 
equation(4) reduces to 

t£(v) + h 2 - v( v j\ nn(v) = o, (8) 



where V(rj) = %- and prime denotes the derivative with respect to ?7.Equation(8) looks like 
the schrodinger equation for a particle having the potential energy V(r/).For high frquency 
waves or for waves such that n 2 >> -^- the general solution to Eq(8) has the form 



Uf 



_L (a n e-^ + (3 n e mr > 



V2 



(9) 



Where a n and /3 n are arbitrary complex numbers. In the opposite limit, n 2 << -^-,the 
general solution is 



H n = A n R + B n R \ 
Here also A n and B n are arbitrary complex numbers. 



*>/ drj 
~R 



(10) 



3 GRAVITON SPECTRA 

It is known that a travelling gravitational wave passing through the barrier V{rf) will always 
be amplified[10].In the quantum treatment of the problem, one says that the initial vacuum 
state of gravitons goes over to a final multi-particle quantum state - the particle creation 
takes place. The final quantum state belongs to a class of the so-called squeezed quantum 
statesfll]. It is clear from Eq(8) that the behavior of u n is determined by the potential 
V(rj) = -^-.The solution(9) is valid outside the potential barrier V(r]) while the solution 
(10) applies to the region where the function V(r]) dominates over n 2 .The values of rj where 
the regimes (9) and (10) interchange are determined by the condition 



n 



\V{rj)\ 



(11) 



A wave with a given "n" enters the barrier region and leaves it in the "turning points" 
defined by eq(ll).We will denote these points by r\i and rjf .The indices "i" and "/" will also 
be used to distinguish the values of the various functions at the initial and final turning 
points. We will match solutions (9) and (10) at rji and rjf for continuity [12,13]. To the left 
of the barrier, that is for r/ < r/i, the solution is taken in the form 
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(12) 



To the right of the barrier, that is for rj>rjf, the solution has the form(9) where the 
numerical values of the co-efficients a n and (3 n follow from the joining conditions. The actual 
values of a n and (5 n are 
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As R{i]) = a(ry)(/> 1 ' 2 (?7) ,one has 
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R' = 6 l l 2 Ha 2 + 
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(15) 



where H = - is used. We write below the complete expression for a n and (3 n . 
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In fact, we have calculated the coefhcients of the Bogoliubov transformation which relates 
the in and out quantum states of gravitons[14]. Let us recall that in the quantum theory the 
quantity |/3 n | 2 defines the spectral energy density of gravitons created from the in vacuum 
state. The contepmorary energy density of gravitons integrated over some frequency interval 
is given by the formula[13,14,15] 

E= I E(v)dv (18) 



The density of states is dN = 
polarization states is given by 



&W 



dv and the corresponding E{v)dv summed over two 



E(y)dv = 87n/ 3 di/|# 
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(19) 



with h=c=l. The remarkable property of Eq(16) and Eq(17) is that they contain only 
the initial "i" and final "/" values of a, H and <^>.Equations(16) and (17) can be further 
simplified for particular numerical values of Hitii, Hfaf and -^-.For instance, in the case of 
(Hidi) 2 ) >> n 2 and (Hfaf) 2 » n 2 ,one obtains[12] 
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On the otherhand,if (i^jOj) 2 ~ 
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,and -f- >> 1 one derives 
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We will consider in more detail the case where the inflationary stage ends with a rapid 
transition to the radiation-dominated stage. The relevant barrier in Eq(8) is described by 
the function V(rj) which first increases with n and then decreases sharply up to zero. Let 
us consider the waves interacting with this barrier. Their wave numbers obey the conditions 
(Hidi) 2 w n 2 ,(Hfdf) 2 » n 2 and we assume also ^- » l.For acceptable inflationary 



models these waves have a present day frequencies in the interval W 8 Hz < v < 10 1& Hz 
[12]. From eq(16) and eq(17) we obtain 
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Using \P n \ in eq(19) with v = #r, we derive 
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The difference between this expresion and the expression in ref[12 ] is the presence of 
contributions due to the scalar field .The kinetic energy of the BD scalar field makes a 
non-zero contribution to E{y) and this comes from the rapid variation of (f> during the 
inflation epoch. When <^>=constant=l, this expression reducess to the expression obtained 
by Grishchuk and Solokhin[12].The extra terms are due to the rapid variation of (f> during 
the inflationary era. 



4 HUBBLE PARAMETER AND MATTER ENERGY DEN- 
SITY 

The entire dependence on n(and hence ,on frequency v) is contained in the factor Hf .In its 
turn, Hf(n) is proportional to the matter energy density at the inflationary stage , namely, 
at those instants of time t n ,when the waves with the corressponding wave numbers "n" 
were entering the under barrier region Hf(n) = ^-p(t n ). The value of H 2 - is propertional 
to the matter energy density at the end of inflation. Some part 5 of this energy((5 ~ 1) has 
been transformed into the electromagnetic radiation. In the course of cosmological expansion 
the energy density of the radiation was decreasing in proportion to (-^-) 4 .The present day 
energy density Ey of the radiation (3°K microwave background) is related to H 2 r according 
to 
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Using the expression for Hf and H ?, neglecting the numerical factors of order unity, and 
restoring p p for the correct dimensionality, we obtain finally 



E(v) = E, 
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This very simple expression allows us to link the energy density E{y) of relic gravitational 
waves measurable today with the time dependent values of the matter density "p" and 
the Hubble parameter "H" attributed to the very early universe .One can say, a little 



loosely,that the function E(v) stores the information on the rate of expansion of the uni- 
verse in that distant past when the relic waves with frequencies v first started to emerge 
from zero-point quantum fluctuations. If the expansion law a(t) were precisely known, 
Eq.(27) would give us a definite spectral dependence of E(v).For example, for the power- 
law scale factors a(t) ~ t a , one has Hit) = at -1 . From the condition Hi(t n )ai(t n ) = n we 

1 — 2 

derive n ~ t®' 1 , that is , t n ~ n^r. Since p(t n ) ~ Hf(t n ) ~ t~ 2 ~ n"- 1 one can obtain 

-2 

E{y) ~ i/a-i , that is the power law frequency dependence for E(v)ln the case of strictly de 
Sitter expansion(i?j =constant) one would obtain the spectrum E(y) = constant. Equation 
(29) gives p and H as functions of frequency v, 

fl» = *£p(y), (31) 

where ^ c is contribution due to the presence of scalar field (j) [square bracket term in 
eq(29)]and, therefore, gives p and H as nonmanifest functions of time, since every v corre- 
sponds to some t n . 

5 CONCLUSION 

We have derived the energy density and spectrum of gravitons in the context of Brans- 
Dicke theory.lt is found that due to the presence of scalar field <fi , the expression for 
graviton energy density is modified [12]. The kinetic energy of <p enters into the expres- 
sion .This happens due to rapid variation of the scalar field 4> during the inflationary 
era . Also the time dependence of the early Hubble parameter and matter energy den- 
sity are derived we find that the frequency dependence of the relic graviton energy den- 
sity spectrum replicates exactly the time dependence of the very early Hubble parameter. 
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